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ABSTRACT 



Context. The presence of titanium oxide (TiO) and vanadium oxide (VO) gas phase species is searched for in the atmosphere of the 
hot Jupiter HD 209458b. 

Aims. We compared a model for the planet's transmitted spectrum to multi-wavelength eclipse-depth measurements (from 3 000 to 
10 000 A), obtained by Sing et al. (2008a) using archived HST-STIS time series spectra. We make use of these observations to search 
for spectral signatures from extra absorbers in the planet atmosphere between 6 000 and 8 000 A. 

Methods. Along with sodium depletion and Rayleigh scattering recently published for this exoplanet atmosphere, an extra absorber 
of uncertain origin, redward of the sodium lines, resides in the atmosphere of the planet. Furthermore, this planet has a stratosphere 
experiencing a thermal inversion caused by the capture of optical stellar flux by absorbers that resides at altitude. Recent models 
have predicted that the presence of TiO and VO in the atmosphere of HD 209458b may be responsible for this temperature inversion. 
Although no specific TiO and VO spectral band head signatures have been identified unambiguously in the observed spectrum, we 
suggest here that the opacities of those molecules are possible candidates to explain the remaining continuous broad band absorption 
observed between 6 200 and 8 000 A. To match reasonably well the data, the abundances of TiO and VO molecules are evaluated from 
ten to one thousand times below solar. This upper limit result is in agreement with expected variations with altitude due to depletion 
effects such as condensation. 
Results. 

Key words, planetary atmospheres - extrasolar planets - HD209458b 



1. Introduction 

The discovery of transiting extrasolar giant planets (EGP) 
has opened the window to the direct detections and charac- 
terization of their atmospheres. Because of the wavelength- 
dependent opacities of absorbing species, measurement of rel- 
ative changes in eclipse depth as a function of wavelength dur- 
ing primary transit has the potential to reveal the presence (or 



(Hubeny et al. 2003 1. In the lower part of the atmosphere (~30 
mbar), the temperature is found to be in the range 1 900-2 400 K 
(Lecavelier et al. 2008b), which corresponds to the M/L/T brown 
dwarf regime, as expected for a hot Jupiter such as HD 209458b 
(|Kirkpatrick"2005]l. 



absence) of specific chemical species (Seager & Sasselov 2000 
IHubbard et al. 2001||Brown 2001| l. In the case ofHD 209458b's 
atmosphere, the transmission spectroscopy method led to the de- 
tection of sodium (Charbonneau et al. 2002). The small atmo- 
spheric sodium signature has made its detection from ground- 
based telescopes difficult. Only recently, Redfield et al. (2008) 
and Snellen et al. (2008) reported the detection of Nai in 
the atmosphere of HD 189733b and HD209458b, respectively 
Absorptions of several percents for H i Lyman-a, O i and C ii 
have been measured in the hydrodynamically escaping upper 
atmosphere ( |Vidal-Madjar et al. 2003| 2004, 2008). More re- 
cently, Rayleigh scattering by H2 molecules has been identified 
(Lecavelier et al. 2008b) and a temperature-pressure (TP) pro- 
file with inversion was derived (Sing et al. 2008b). This temper- 
ature inversion leads to high temperature at both low and high 
pressure. Note that this temperature bifurcation was very well 
predicted by atmospheric models of strongly irradiated planets 



In the cool atmosphere of sub-stellar objects like brown 
dwarfs and close-in EGPs, the strongest absorption fea- 
tures are expected to be those of alkali metals like sodium 
and potassium ([Seager & Sasselov 2000 Sudarsky et al. 2000 



IBurrows et al. 2000 i This is due to the condensation and the rain 
out of elements in the atmosphere which clears the atmosphere 
of most of the metals and keeps the less refractory alkaU met- 
als (BuiT ows & Sharp 1999, Lodders 19991 IBurrows et al. 20021 
Sudarsky et al. 2003| l. Thus, in contrary to the spectrum of 
HD 189733b (Pont et al. 2007; Lecavelier des Etangs et al. 
2008a), the spectrum of HD 209458b should be dominated by 
absorption from sodium Nai and potassium Ki. 

However, depending on the effective temperature, a large 
number of diatomic and polyatomic molecules are predicted 
to be present according to various models of brown dwarf 

ILodders 19991 



( [Burrows & Sharp 1999 
ILodders 20021 



Hubeny et al. 2003 



Send offprint requests to: desert@iap.fr 



and hot Jupiter 

All ard et al. 200n 

Bur rows et al. 2006] l. Among those molecules and at a 
temperature above 1 800 K, titanium oxide (TiO) and vanadium 
oxide (VO) in gas phase equilibrium are most probably present 
with a high abundance in strongly irradiated planet atmospheres 
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(Seager et al. 1998; [Hubeny et al.^003| [Fortney et al.lMi} . 
Furthermore, the low albedo measurements of HD 209458b 
dRowe et al. 2 006 ) rule out most of the absorbents, but TiO 
and VO. Recently, the interpretation of visible observations 
obtained during the transit of this planet (Sing et al. 2008a), as 
well as the interpretation of near infrared observations taken 
during its secondary eclipse (Knutson et al. 2008) require that 
HD 209458b experience an inversion in the atmospheric P-T 
profile and a have stratosphere (Sing et al. 2008b; Burrows et 
al. 2007c; Burrows et al. 2008). Such an inversion could be due 
to the absorption of the incident visible light by TiO and VO 
molecules (Fortney et al. 2007). Ongoing models (Burrows et 
al. 2008) claim that this inversion must be caused by the capture 
of the incident optical stellar flux by a stratospheric absorber 
that could potentially be TiO and VO. Similarly other studies 
( [Hubeny et al. 2003 Fortney et al. 20071 have more particularly 
investigated the effects of TiO and VO which absorb much of 
the incoming stellar radiation high in the atmospheres of the 
hottest CEGPs. 

In this work, we use STIS spectra obtained during plane- 
tary transit at two spectral resolutions (low and medium) to 
search for the direct evidence of the presence of TiO and VO 
in the planet's atmosphere. Both datasets are combined to ex- 
tend the measurements over the entire optical regime to quan- 
tify other possible absorbers appearing in the transmission spec- 
trum. Rayleigh scattering and sodium absorption have been pro- 
posed to explain the spectral features observed between 3 000 A 
and 6 000 A (Lecavelier des Etangs et al. 2008b; Sing et al. 
2008a; 2008b) yet, additional absorptions remain unattributed 
from 6 100 A up to 8 000 A. These absorptions cannot be due 
to optically thick high altitude clouds, which would otherwise 
mask the detected signature of Rayleigh scattering (Lecavelier 
des Etangs et al. 2008b). Since TiO and VO molecules are ex- 
pected to be present in high abundance for the range of tempera- 
tures and pressures for strongly irradiated hot Jupiter, their opac- 
ities could largely contribute to the observed absorption in this 
wavelength domain (Seager et al.l998; [Burrows & Sharp 1999{ 
Sharp & Burrows 2007| l. 



After a brief description of the observation and interpretation 
status (Sect. |2]i, we present, in Sect. [3l the model together with 
the method used to calculate the atomic and molecular opaci- 
ties. Finally, we estimate the contribution of TiO and VO to the 
observed transmission spectrum for several abundance scenar- 
ios (Sect, m and discuss the consequences on the atmospheric 
chemistry and aeronomy of HD209458b (Sect.jSj. 



2. Observations 

The HST-STIS G750L, and G430L low resolution grating ob- 
servations of HD 209458b analyzed here are also detailed in 
[Knutson et al. (2007a)| [BaUester et al. (2007)[ and Sing et al. 
(2008a). For both flie G750L, and G430L gratings, two vis- 
its were done for each grating, of five consecutive orbits each. 
Each visit span one completed transit of HD 209458b in 
2003. Together they cover the combined range 2900-10300 A, 
with some overlap around 5 300 A with a resolving power 
R = 500. While the broad-band spectrophotometric analysis of 
[Knutson et al. (2007a)[ was used to identify water absorption fea- 
tures in the atmosphere of HD 209458b at wavelength greater 
than 9 000 A JBarman 20071) . here we analyze the spectrum at 
full resolution. 

We directly use HD 209458b absorption depths (AD, see 
definition in Sect. [33) between 4000 A and 8 000 A pro- 



vided by [Sing et al. (2008a)| The low resolution transit spec- 
trum ratio, limb darkening corrected, is plotted in Figure [T] 
The spectrum below 8 000 A is considered here where the ab- 
sorption due to water molecules is negligible. Three promi- 
nent broadband absorption features are revealed. The near ul- 
traviolet (NUV) absorption, at wavelength A < 5 000 A, was 
first reported by [Ballester et al. (2007)[ and explained by the 
absorption of a hot hydrogen layer within the atmosphere. 
[Lecaveher Des Etangs et al. (2008b) propose an alternative ex- 
planation invoking the Rayleigh scattering by H2 molecules. 
Within the same datasets. Sing et al. (2008b) found that the Na 
spectral line profile is characterized by a wide absorption with 
a sharp transition to a narrow absorption profile at higher alti- 
tudes values. This sharp transition is interpreted by condensa- 
tion or ionization which deplete Na atoms in the upper atmo- 
sphere. Using a global fit to these data, from 3 000 A to 6 200 A, 
[Sing et al. (200 8b) determine the average pressure-temperature 
profile (P-T) at the planetary terminator Two types of P-T pro- 
file were derived which both can lead to the depletion of Nai 
atoms. One profile implies the condensation of Na i and the other 
the ionization of Na i. We use this former one (Fig. [2]). Both pro- 
files experience hot temperatures, above 1 800 K at 50 mbar, due 
to Rayleigh scattering (Lecavelier et al. 2008b). At such a tem- 
perature, TiO and VO are in gas phase. 

A third broad absorption feature in the range 6 200 - 8 000 A 
appears on the AD spectrum and is still unexplained. This ab- 
sorption is composed of two remarkable features visible in 
Fig. [I] 

1. A strong and broad spectral absorption feature centered at 
about 6250 A which is slightly smaller than for the Nai 
(AD - 0.0148% over 100 A). Other spectral signatures ap- 
pearing at higher wavelengths are weaker than the 6 250 A 
one. 

2. A broad and continuous absorption level above 7 000 A 
systematically higher than the lowest observational AD 
(0.0146% over 1500 A). This broad continuous flat absorp- 
tion excess is weaker than the spectral feature at 6 250A but 
is equally as important since it covers a large spectral do- 
main, and has a much larger absorption level than predicted 
by extrapolating Rayleigh scattering and Nai or Ki broad 
wings. 

We overplotted in Fig. [T] the radii presented by Knutson et 
al. (2008) that were derived using spectrophotometric bandpass 
across roughly thousand angstrom bins. These radii are found to 
be slightly different than those presented in this paper, especially 
in the region of interest of this study (over 6 000 A), though 
the differences are either close to or within the error bars of the 
two studies. This differences could be explained by the supple- 
mentary corrections in the data treatement employed by Sing et 
al. (2008) who corrected the transmission spectrum for telluric 
contamination and various systematics, quantified the level of 
systematic red noise to take it into account in the computed er- 
rors and finally applied a spectral limb darkening correction be- 
fore binning. Note that the systematic errors have a weak wave- 
length dependance, affecting the absorption depth values over 
large wavelength regions in a similar manner, largely preserving 
the wavelength structure of narrow bands such as the 6 250 A 
feature. 

Looking at the low resolution spectrum, we find the mini- 
mum AD of 0.01453% around 5 000 A, which indicates a base 
pressure of 33 mbar in the P-T profile(Lecavelier des Etangs 
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et al. 2008b). Conversely, the maximum AD of 0.0149% is ob- 
tained for the central pixel of the Na i line which corresponds to 
a pressure of 1 mb (Sing et al. 2008b). 

Below we present possible identifications of absorbers that 
can simultaneously explain the spectral features between 6 200 
and 7 000 A and the broad relatively flat absorption excess from 
7 000 to 8 000 A. 

3. Model description 

In order to interpret the transit observations, we developed a 
model that calculates the absorptions of atomic and molecular 
lines in the atmosphere of a transiting planet. 

Transit spectroscopy probes the transition region between 
the day and night sides, called limb. We used the geometry of 
a transiting system (Br own 20011 1 and a model originally de- 
veloped by Ehrenreich et al. (2006) and adapted for this work. 
Along a single cord, stellar photons cross several levels of the 
spherically stratified atmosphere of the planetary limb. 

The total opacity, t,;, along a cord (parallel to the line of 
sight) as a function of its impact parameter b and the altitude h, 
is given as: 



TAAb) = 2 



f 

Jo 



o-A,ini(h)dl, 



(1) 



where cr^j is the cross section coefficient for the species / at the 
wavelength A, and «,(/;) is the density of the species / at an alti- 
tude h in the atmosphere (see Fig. 1 in Ehrenreich et al. 2006). 

3.1. Atomic and molecular line opacities 

textbfin the case of alkali atoms, we calculate the opacity us- 
ing the oscillator strength, the radiative lifetime and the colli- 
sional broadening taken from Morton et al. (1991) and Iro et al. 
(2005) which used Burrows et al. (2000) for determining the line 
shapes. Further progresses have been achieved in the theory of 
these profiles (Burrows et al. 2003; AUard et al. 2005). The dif- 
ference between these new line shapes and the ones used in this 
study resides mainly in the far wings. However, the contribution 
of other absorbents dominates the spectrum in the far wings of 
the alkali lines. Thus, we found that the calculations of Morton 
et al. (1991) and Iro et al. (2005) were precise enough for our 
study. 

In the case of molecules, we calculate the line strengths 
for each line of each species using the available data. We 
derive the strength at local thermal equilibrium (LTE), noted 
S, of the molecular spectral line from the general equation 
( [Sharp & Burrows 2007| l: 



nieC Q{T) 



,-hc(Fj-F,)/kT] 



(2) 



where gi and fij are the statistical weight of the / energy level 
and the oscillator strength for a transition from that level to a 
higher level j, respectively. The excitation energies F/, and Fj 
are the term values in cm"' of the and /'' levels participat- 
ing in the transition. In addition, Q{T) is the partition function 
of the species at some temperature T. We compute the partition 
function as described by |Sauval & Tatum (1984)[ considering a 
fourth degree polynomial function of log(^Y'^). The molecular 
lines are broadened using a Voigt profile. At pressure P, the col- 
lisional broadening is calculated using (Sharp & Burrows 2007): 



with Wo = 0.1 cm"' bar"' 

Finally, we obtain the total monochromatic opacity of the 
atmosphere by summing the individual contributions for each 
atom and molecule weighted by their respective abundances at 
each altitude level. 



3.2. Grid of opacities 

Computing the opacities directly at each altitude with the cor- 
responding temperature and pressure would require intense cal- 
culation. Instead, we precalculate a grid of opacities for each 
molecule. Our grid is composed of 42 cells with a range of tem- 
perature from 100 to 5 000 K and a range of pressure from 1 to 
10"^ bar. Opacities have been computed from 2900 to 10000 A 
with a resolution of 0.005 A which corresponds to 5 million 
elements per cell. For a given P-T profile we then calculate the 
corresponding opacities by applying a linear interpolation on the 
grid cells. 



3.3. Absorption Depth (AD) 

The output of the model is a spectrum ratio defined by Brown 
(2001) as the ratio of the flux received during the transit with the 
flux received when the planet is not occulting the star of radius 
R-t. Otherwise indicated by the absorption depth (AD), that is 
the occultation of the surface of the star by the surface of the 
planet at each wavelength. The surface of the planet is twofold: 
the optically thick disk of radius Rp and the wavelength- 
dependent occultation by the atmosphere that surrounds the 
planetary disk expressed as an equivalent surface. We used 
Rp = 1.32 Rjup andRj, = 1.1 25 Rq from Knutson e t al. (2007a)| 



Lecavelier Des Etangs et al. (2008b)^ show that the observed AD 



is well approximated by: 



2H ax 
AD, =AA)(1 + — In— ) 

Rp cr^o 



(4) 



where ADo is AD sS. X - Aq and H the scale height. Thus the 
observed mean AD over a wavelength range of a given spectral 
element (typically 50 A), is proportional to the mean of the log- 
arithm of the cross section. Therefore, we calculate the eff'ective 
cross section in a given spectral element by averaging the loga- 
rithm of the cross section calculated at a much larger resolution: 



(5) 



Av, 



Lor en t: 



= WoP, 



(3) 



4. Analysis 

Spectral signatures, line profiles and abundances vary de- 
pending on the temperature and pressure vertical profile in 
the atmosphere. Here we use the vertical P-T profile pre- 
sented in Figure |2] and two sodium mixing ratios from 
Sing et al. (2008b) as well as the Rayleigh scattering derived by 
Lecaveh^ erPes Etangs et al. (2008b)| The AD of 0.01453%, cor- 
responding to a pressure of 33 mbar, is taken as the reference at 
1.32 Jupiter radius (Knutson et al. 2007). Using those hypoth- 
esis, we calculate the absorption spectrum with only the abun- 
dances as free parameter We vary the abundances to find the 
best fit to the data using ax^ minimization between the data and 
the model. 



4 



J.-M. Desert et al.: TiO and VO broad band absorption features in the atmosphere of HD 209458b 



Table 1. List of the most probable species with their corre- 
sponding solar abundances which have significant absorptions 
in the 6000-8 000 A wavelength domain (from Lodders 2002 



and Sharp & Burrows 2007 1. This list is ordered by the product 
of an average of the cross section between 4 000 and 8 000 A 
and its corresponding solar abundance. This ordered list gives a 
rough estimate of the plausible detectability for each atomic and 
molecular species. 



Element 


Xa 


< cr,, > (cm^) 




Nai 


3.e-6 


Le-19 


l.e-25 


Ki 


2.e-7 


Le-20 


l.e-27 


Li I 


Le-9 


Le-19 


l.e-28 


SiO 


Le-5 


l.e-17 


l.e-22 


TiO 


Le-7 


l.e-16 


l.e-23 


VO 


Le-8 


l.e-16 


l.e-24 


MgH 


Le-9 


l.e-17 


l.e-26 


H2O 


Le-3 


l.e-24 


l.e-27 


FeH 


Le-9 


l.e-19 


l.e-28 


CaH 


l.e-11 


l.e-17 


l.e-28 


CrH 


l.e-10 


l.e-19 


l.e-29 



4.1. Atomic line: 

The most important absorptions by atomic lines are those of al- 
kali metals, Nai, Ki and Lii (see Table [TJ. This table provides 
the lines with the strongest abundance weighted cross-section in 
the considered wavelength range. None of this atomic lines, nor 
Ha, can explain these observed spectral feature at 6 250 A. 

In the low resolution transit spectrum, we do not detect sig- 
nificant absorption in the K i line at 7 698 A nor in the Li i doublet 
at 6 708 A. These non detections exclude the possibility of so- 
lar abundance over the whole atmosphere for these two species 
(see Fig.[Tll- Assuming a constant mixing ratio, we constrain the 
Ki and Lii abundances to be lower than 2 x 10"^ and 2 x 10 ' 
solar, respectively (Icr). Of course, these upper limits does not 
exclude larger abundances below the altitude corresponding to 
the observed absorption depth of ~0.0145%. 

We conclude that the absorption depth measured between 
6 200 and 8 000 A cannot be explained by the strong absorption 
lines from abundant atomic species. The line profile of Nai does 
not affect the results on TiO/VO abundances since the remaining 
broad band absorption begin at 200 A away from the center of 
the line on the red side. In the following, we consider possible 
absorptions by molecules. 

4.2. Molecular line: 

A list of plausible molecular species absorbing in the 6 000 A 
region is presented in Table [1] 

Although the molecules CH4, NH3 and CO are the most im- 
portant source of opacities in the infrared, their effects are negli- 
gible in the visible wavelength domain; they are not considered 
in the following analysis. 

Using the abundance weighted mean cross-section list, we 
found that most plausible species are SiO, TiO and VO respec- 
tively. Silicium oxide (SiO) absorption mainly occurs below 
2 000 A, outside our wavelength domain. As a result, TiO and 
VO are predicted to be the most abundant and the most important 
sources of opacity in the present domain ( Fortney et al. 2006b , 
ILodders 20021) . However, since both species condense, their 
abundances decrease below 1500K. 



Amongst all the other molecular absorbers that can be used 
to explain the observation in this wavelength range, the hydride 
metals CrH, FeH, MgH, and CaH were considered as poten- 
tial candidates to explain the shape of the spectrum. However, 
the abundance of those elements drops with decreasing tem- 
perature below 2 000K, due to the formation of condensates 
(Lodders & Fegley 2006 Kirkpatrick 20051). Furthermore, their 
absorption domain does not match the one considered here, like 
for H2O. 

Among plausible molecules in the atmosphere of 
HD 209458b, TiO and VO have the largest cross-section 
in the visible domain. We explore the extent to which var- 
ious abundances of TiO and VO can be used to reproduce 
both the medium and low resolution spectra. We derived the 
LTE strength S (Eq. |2]) of the molecular spectral lines from 
|Plez (1998)1 These data Ust about 3 million lines from 4000 
to nearly 30,000 cm covering the whole visible part of the 
spectrum for the transitions of the most abundant isotopes '^^TiO, 
together with ^^TiO, "^^TiO, '^'^TiO, and ^"TiO, for nine different 
electronic band systems. Using the grid previously described in 
Sect. 13.21 we compute theoretical AD including TiO and VO 
absorption at each wavelength (see Fig. 13). At a temperature 
of 2 000 K and a pressure of 0. 1 bar, the lines are significantly 
broadened, thus the rapid fluctuations in cross-section over short 
wavelength intervals are suppressed revealing the main band 
features. At significantly lower pressures, the broadening of the 
lines is much smaller, the cross-section shows rapid variations 
as a function of wavelength and the broad main features do not 
show up so clearly. 

In the following sections, we explore different scenarios for 
the atmospheric abundance of TiO/VO in order to adjust our 
model to the observations. We fit the data with an atmospheric 
model including Na I, Rayleigh scattering, and TiO and/or VO 
in different abundances. In a first step, we will consider a unique 
abundance of TiO/VO over the whole atmosphere. In a second 
step, we will assume two levels of abundance of TiO/VO. We fi- 
nally discuss the fits obtained using those different abundances. 

4.3. TiO : 

As a first step, we consider only TiO and constrain the abun- 
dance by fitting the observational low resolution spectrum ratio. 
Starting with the simplest assumption, we consider a solar abun- 
dance over the whole atmosphere of the planet. 

Our model gives an AD average of 0.155% over the 
full bandpass 4000-8000 A (see upper dashed-dotted line in 
Fig. Eh)- Such a value is significantly above what is observed. 
Therefore, TiO abundance must be much lower than solar above 
altitudes corresponding to the observed AD level. 

Assuming a constant mixing ratio, the best fit of the low res- 
olution data set over the region 4 000-5 700 A gives a TiO abun- 
dance of (8. ± 0.6) X 10""^ solar. With this TiO abundance, the;^^^ 
is reduced from 1754 to 1207 for n - 942 degrees of freedom 
(using non rebinned spectra). However, this abundance cannot 
reproduce the 6 250 A spectral feature and the flat absorption of 
the observations between 6 200-8 000 A. We conclude that other 
absorbers must be considered to explain this part of the AD spec- 
tral profile. 

4.4. TiO and VO: 

Since VO is the most important source of opacity after TiO (see 
Table [U, we add this molecule to our model. The chemistry of 
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the vanadium is quite similar to the one of titanium, where the 
monoxide is in the gas phase at a temperature higher than the 
temperature of condensation. 

In a first step, we consider a solar abundance over the whole 
atmosphere. As for TiO, VO with solar abundance cannot re- 
produce the observations, especially around 5 000 A. The abun- 
dance of VO must be lower. Assuming constant mixing ratios for 
both molecules, the best fit gives for TiO (6 + 0.6) x lO"'* and for 
VO (3 + 0.5) X 10"^ as shown in Fig.|4j5. We obtain a satisfac- 
tory fit with ax^ of 817, or ax^/n ^ 0.86 in the low resolution 
data set. The continuous flat part of the AD curve can be well 
reproduce by absorption of VO molecules. The region around 

5 000 A already well reproduced with Rayleigh scattering and 
with sodium lines absorptions, is not affected by the presence of 
TiO and VO molecules. However, the spectral features around 

6 250 A remains unexplained assuming constant mixing ratios 
for both species. Note that sligthly lower abundances are nec- 
essary for TiO and VO if we consider the Knutson et al.(2007) 
results for the radii over 6 OOOA. 

4.5. TiO, VO, and condensation: 

With a large temperature gradient, TiO and VO abundances 
are expected to vary along the vertical atmospheric profile. 
Furthermore, the P-T profile derived by |Sing et al. (2008b)| 
crosses the TiO and VO condensation curve (See Fig. |2]l. To 
check if signature of condensation can be found in the present 
data set, we assume in this section that TiO/VO molecules can 
be depleted above a given altitude. Indeed, TiO and VO are ex- 
pected to condense into many Ti and V-bearing compounds as 
described in details by Lodders (2002). To implement this as- 
sumption, we introduce a maximum altitude, corresponding to 
an absorption depth of condensation (ADcond), above which the 
absorption by the corresponding molecules in gas phase van- 
ishes. 

In the model, we derive the altitude of condensation from 
the intersection between the P-T profile and the condensation 
curve. The condensation curves of TiO and VO depend mainly 
on the temperature (see Fig. [2| and [Sharp & Burrows 2007| ). The 
TiO condensation point occurs at a temperature of I 850 K and 
a pressure of 0.02 bar (see Fig.|2]l. This altitude corresponds to 
an AD of condensation ADcond - 0.0146%. In the case of VO, 
the condensation point is at 1 600 K and 0.01 bar, leading to 
ADcond = 0.0147% (see condensation limits in Fig. [4]:). 

Using the P-T profile presented in Fig. |2l TiO and VO 
ADcond levels are lower than the 6 250 A spectral feature AD 
peak. However, both ADcond levels are found to be above the 
observed average AD, above the continuous flat part between 
6 500 A - 7 500 A, and significantly above the lowest AD - 
0.0144% around the left wing of the Nai at 5 000 A. Thus, 
these molecules cannot be in solar abundances below the alti- 
tude of condensation, otherwise the resulting AD curve would 
be at the level of the condensation limits. We fit the AD curve 
between 4 000-8 000 A and obtain an abundance for TiO and 
VO below their altitudes of condensation of (8 ± 0.5) x lO"** and 
(6 + 0.4) X 10"^ solar, respectively. We obtain a fit with a similar 

as before (Sect. |4~4l ) with one level of abundance for TiO/VO 
over the whole atmosphere. 

Although the "cold-trap" can explain the lowest observa- 
tional ADs, our assumption of null abundances above this al- 
titude cannot describe the observed 6 250 A spectral signature. 
Thus TiO/VO molecules should remain in the upper part of the 
atmosphere if the feature is due to these species. 



Table 2. Bests xi fits performed for various model and the de- 
rived TiO and VO abundances. 



Model 


X2 


TiO(xl0-* 


VO (xlO--) 


Only TiO 


1207 


(8. ±0.6) 





TiO and VO 


817 


(6. ± 0.6) 


(3. ± 0.5) 


TiO and VO with condensation 


816 


(8. ±0.5) 


(6. ± 0.4) 


2 levels of TiO and VO 


802 


(1. ±0.2) 


(1. ±0.2) 



4.6. TiO and VO in two separated iayers: 

In this section, we explore the possibility of the presence of TiO 
and VO molecules in the upper part of the atmosphere above 
the condensations altitudes. In the lower part of the atmosphere, 
we keep the abundances found in Sect. 14.51 Then a level of null 
abundance is imposed, corresponding to the condensation level. 
Finally another abundance level is chosen for the upper part of 
the atmosphere, where the temperature rise up and cross again 
the condensation curves of TiO/VO. In the upper part of the at- 
mosphere, where the pressure is lower, the TiO lines are less 
pressure broadened, numerous narrow peaks appear in the TiO 
opacity curve, in particular, a strong peak appears at 6 200 A 
(Fig. Si). 

The best fit gives a TiO and VO respectively (1 ± 0.2) x 10""* 
and(l+0.2)xl0"~ solar abundances in the upper layer, above the 
altitudes of condensation. However the x^ does not significantly 
decrease (802 for 944 degrees of freedom). The spectral features 
of the AD curve is now better fitted with TiO, at least to first 
order, whereas VO molecules reproduce mainly the continuous 
flat part. In that way. We can reproduce more successfully the 
spectral feature, the continuous flat part and match the Na i blue 
wing simultaneously. Note that the flat part of the spectrum ratio 
cannot be reproduced if we choose a null abundance of TiO/VO 
in the lower part of the atmosphere. 

5. Discussion 

We have highlighted that the HD 209458b atmosphere is op- 
tically thick at low pressures and this requires new absorbers 
between 6 000-8 000 A. The observation of the Na i in the atmo- 
sphere exclude the presence of clouds at pressure under 30 mbar 
in the observed limb of the planet. We have shown that Li i and 
Ki absorption cannot reproduce the broad band absorptions ob- 
served. Nonetheless, we derive upper limits for their abundances 
at 2 X 10-^ solar for K i and 2 x 10"' for Li i. 

The Barman (2007) analysis of the Knuston et al. (2007) re- 
sults and Sing et al. (2008b) analysis both agree that remaining 
absorbents are required to explain the continuous broad band ab- 
sorption observed between 6 200 and 8 000 A, since the radii 
derived in both reductions are significantly above the model 
with no TiO nor VO molecules (see Fig. As a consequence, 
the presence of remaining absorbents, possibly TiO and VO, is 
mandatory to explain the observations. This result does not de- 
pend on the P-T profile adopted. The abundances of TiO and VO 
are only poorly constrained by any P-T profile. We decided to 
use the P-T profile derived by Sing et al. (2008) as an example. 
In that case, the lower part of the P-T profile, i.e below the con- 
densation curves, is the most critical part for the determination 
of the TiO and VO abundances. The abundances of TiO and VO 
molecules are mainly constrained by the level of the observed 
AD curve between 6 500 and 8 000 A and by the pressure of the 
lower point. Indeed, for this lower point, the TiO and VO lines 
are mostly pressure broadened. Thus, the abundances of these 
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species are mainly sensitive to the pressure and not to the tem- 
perature in the lower part of this profile. The P-T profile above 
the condensation curves constrains mostly the spectral features, 
especially the 6 250 A one, as seen for the case two level of 
TiO and VO. However, we do not detect any TiO/VO narrow 
spectrale features, thus any P-T profile above the condensation 
curve can be adopted. With the P-T profile used in this study, we 
have shown that no clear solutions emerge wether condensation 
is taken into account, or if two levels of TiO and VO molecules 
are considered (see Table |2]i. Thus, other profiles could led to 
the same conclusion for the presence of TiO and VO and/or of 
potential additional absorbers. 

Although no typical TiO and VO spectral signatures have 
been identified unambiguously in the observed spectrum, we 
suggest that the opacities of those molecules are the best can- 
didates to explain the remaining continuous broad band absorp- 
tion observed between 6 200 and 8 000 A. Theoretical absorp- 
tion with models including TiO/VO with abundances below so- 
lar were evaluated to match the data reasonably well. Using the 



P-T profile from Sing et al. (2008b)] we derived upper limits for 
the TiO and VO abundances. The model without TiO or with 
TiO but no VO molecules give the worst ;i'2 solutions. We found 
that the abundance of TiO should be around 10 to 10"^ so- 
lar, and the abundance VO around 10"^ to 10"^ solar. The fits 
become marginally worse on the blue side of the sodium wing 
when adding TiO and VO (see Fig.|4]i. This is because the P-T 
profile used in our work was obtained by fitting the observation 
on the blue side of the sodium line with a model which contains 
only Rayleigh scattering and sodium opactity (Sing et al. 2008b) 
which helped limit the number of free parameters in the fit, such 
that a convergent solution could be found easilly. This model fits 
precisely the observed data for this part of the spectrum (see the 
dotted line in Figs. [T]and |4]i. The addition of the TiO and VO 
opacities in this model, contributes to slightly decrease the qual- 
ity of the fit on the blue side of the sodium line (see continous 
line in Fig.|4]i. This is, however, of marginal consequence on the 
TiO and VO evaluations in our study since these are mostly con- 
trolled by the main part of the observations in the redder range 
of the spectrum. Considering the quality of the observations this 
does not affect the TiO and VO evaluations. Hence, the use of the 
P-T profile with only Rayleigh and sodium calculated by Sing 
et al.2008 is sufficient to derive reasonably good estimates of the 
TiO and VO abundances. 

Two layers of TiO/VO in the atmosphere are better suited 
to explain the 6 250 A spectral feature, however the is not 
significantly improved (Table|2|i. Note that, the 6250 A spectral 
feature is seen in the two independent spectra obtained at both 
low resolution, as mentioned here, and at medium resolution as 



observed by Sing et al. (2008a) in the G750M STIS spectra 



The tholins, polyacetylenes, or various non-equilibrium 
compounds at high-altitude could also be responsible for the 
temperature inversion. However, the study of these molecules 
requires a full non-equilibrium chemical model which is not 
the purpose of that work. TiO and VO molecules are very 
good high-altitude absorber candidates (Hubeny et al. 2003 



Fortney et al. 2007). Fortney et al. (2008) highlighted the impor- 
tance of gaseous TiO and VO opacity in their model of highly 
irradiated close-in giant planets. They define two classes of ir- 
radiated atmospheres. Those which are warm enough to have 
a strong opacity due to TiO and VO gases ("pM Class" plan- 
ets), and those that are cooler ("pL Class" planets) dominated 
by Nai and Ki. Our possible detection of TiO/VO and Nai 
in the hot atmosphere of HD 209458b confirm that this planet 



is located in the transition region between the two classes de- 
fined by these authors. As observed here, due to the presence 
of TiO/VO, HD 209458b's atmosphere absorbs incident energy 
between 6 500 A and 8 000 A. Consequently, this planet has a 
hot stratosphere (around 2 500 K) with a temperature inversion 
(IBurrows et al. 2007cl l. As another consequence, because of the 
thermal emission of the energy trapped by the TiO/VO absorp- 
tion, the planet appears very bright during the mid infrared sec- 
ondary eclipse dKnutson et al. 2008] l. 

The first titanium condensates appearing at high 
temperature and high pressure are Ti^Os and CaTiOa 
(Burrows & Sharp 1999 ILodders 2002l l. The condensation 



of vanadium starts at lower temperature ( 1 600 K) than for 
titanium (1 800 K). The vanadium condenses into solid VO 
and then into V2O3. Thus, depletion of TiO should start at 
lower altitude than for VO. Below 1 600K, TiO, VO and major 
refractory elements are absent, leaving monoatomic Na i and K i 
to dominate the spectrum. Nevertheless, sodium and potassium 
chlorides become increasingly abundant with decreasing tem- 
perature, especially for KCl which is the dominant K-bearing 
compound. Condensed potassium depletes the atmosphere of 
atomic Ki, as seen with Nai, leading to reduced signatures. 
This could explain why Nai is seen in abundance in the low 
resolution data, but not Ki. In the lower warmer atmosphere, 
where wide atomic K line wings would be observable even at 
low resolution, TiO and VO likely make up the surrounding 
continuum further masking the signature. 

An alternative P-T profile with Na ionization has been also 
proposed (Sing et al. 2008b). A model with Na ionization leaves 
a wide range of temperatures possible in the middle atmo- 
sphere. Within the framework of equilibrium chemistry, this pro- 
file presents a lower temperature gradient (Fortney et al. 2003). 
The minimum temperature is below the TiO/VO condensation 
curves, and above the sodium condensation curve. We did not 
enter into a detailed study of TiO/VO abundances with these P- 
T profile; however, we note that these profiles would have the 
effect to increase the abundance of TiO and VO molecules in 
the middle atmosphere. See Sing et al. (2008b) for more details 
about ionization. 

If present in the upper part of the atmosphere, 
TiO/VO molecules are not fully depleted by condensa- 
tion. Thermochemical equilibrium calculations with rainout 
ILodders 20021 



(Burrows & Sharp 1999 



Hubeny et al. 2003 1 



have shown that TiO and VO can exist at high-T/low-f points 
such as in the upper part of our P-T profile. The presence of TiO 
and VO at high-T/high-P points leads to a situation in which 
there is two levels of TiO/VO. This has been akeady proposed 
by Hubeny et al. (2003) to explain the presence of a temperature 
inversion in strongly irradiated planets atmosphere's. Since the 
P-T profile used here crosses twice the Ti and V condensation 
curves, two levels containing TiO/VO molecules in gas phase 
could be separated by a middle level that is free of those 
molecules, as proposed by Hubeny et al. 2003. Nonetheless, a 
cold trap region is usually expected to deplete the upper low-P 
region. Flushing out the upper atmosphere of TiO/VO would 
rule out 6 250 A feature as being due to TiO. However, for now, 
it has not been proven that cold-trap effect would deplete those 
molecules, especially for the hydrodynamically escaping type 

2004, 



of upper atmosphere ( Vidal-Madjar et al. 2003 



2008; 



Lecavelier 2004; Lecavelier et al. 2007; IEhrenreich et al. 20081 1. 

New observations with higher signal-to-noise ratio and bet- 
ter resolution, together with improved chemical models, are re- 
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quired to address TiO and VO spectral absorption features and 
detect specific band heads. 
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Fig. 1. The low resolution STIS measurements of the planetary transit absorption depth (AD) corrected from limb-darkenning ef- 
fects and binned by 60 pixels (histogram). The observed Icr errorbar is plotted above the spectrum. The radii derived by Knutson 
et al.(2007) with the corresponding error bars and spectrale bins are plotted (squares). The dashed line corresponds to the best 
fit model assuming Rayleigh scattering and sodium absorption with the physical T-P profile plotted in Fig. |2] The difference be- 
tween the dashed curve and the observed AD spectrum for wavelength over 6 000 A indicates that remaning absorbents should be 
presents in this spectral domain. Overplotted in continuous line, is the same model, with atomic lithium and potassium with solar 
abundances. This model cannot reproduce the observations above 6 200 A. 
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Fig. 2. The atmospheric Temperature-Pressure profile {P-T profile) with error bars derived by Sing et al. (2008b) used in this study. 



The dashed lines correspond to the condensation curves for titanium, vanadium and sodium. The three points of this T-P profile 
are derived from the fit of the observed absorption depth curve. The hot point at a pressure of 0.05 Bar is imposed by the Rayleigh 
scattering (See Lecavelier et al. 2008). 
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Fig. 3. The logiQ of the monochromatic cross-section cr (cm^) as a function of wavelength for the vibration-rotation transitions of 
TiO and VO. The contribution due to different isotopes is included. TiO has a strong absorption feature shortward of 7 500 A, and 
has a strong peak near of 6 200 A as observed in Fig.[T] 
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Fig. 4. Same as Fig. [T] with fits obtain with various models (continuous thick hne). The observed Icr error bars are plotted around 
7 800 A (see also Fig. 10 of Sing et al. 2008). The dashed line shows the fit obtained using only the model with Rayleigh scattering 
and sodium absorption with the T-P profile (Fig.|2]i. A remaining broad band absorption appears in between 6 200 to 8 000 A which 
could be attributed to TiO and VO opacities, a: Model with a constant mixing ratio of TiO. A constant solar abundance of TiO 
is excluded (upper dashed-dotted line). The best fit gives an abundance of (8 + 0.6) x 10"^ solar (red solid line), b: Model with a 
constant mixing ratio of TiO and VO. Solar abundance for VO is also excluded. The red solid line is the best fit which includes TiO 
and VO. c: TiO, VO, and condensation. The altitudes of condensation for titanium and for vanadium are both plotted (dash dotted 
lines) together with the best fit (red solid line), d: Best fit for two distinct levels containing TiO and VO (See Sect. 14. 6t . No clear 
solution emerge in between the three last models (red solid line). 



